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1. INTRODUCTION 

AMB is one of the magnetic suspension methods which enables rotor shaft in rotary machines to be 
lifted off without mechanical support by actively controlling the electromagnet [1]. Absence of mechanical 
contact results in friction reduction. Thus, speed or acceleration of supported components can be increased 
significantly when comparing with conventional bearing in use. Authors in [2] has listed many advantages of 
AMB such as the absence of lubrication and contaminating wear, high speed rotation, low bearing losses, etc. 
These benefits allow AMB to be integrated into electrical motor which can rotate upto 200000 rpm [3] and 
be applicable as ultra-high speed spindle in machine tools [4, 5]. AMB can also be adopted in vacuum and 
cleanroom systems [6] or equipment with harsh working condition like turbo machinery [7, 8]. Design and 
different structures of magnetic bearings have been presented in [2, 9]. It can be seen that the system is highly 
non-linear, and can become very complex. Thus the challenge lies in developing control scheme of bearing 
so that it assures high performance features, especially nanometer accuracy [10] since the gap between rotor 
shaft and bearings can be extremely tiny. Poor design of controller may result in rotor unbalances and internal 
damping which in turns create vibration in machines, crack in motor shaft and failure at the end [11, 12]. The 
classical PID algorithm have been widely used to control AMB system due to its simplicity, adaptability and 
maturity [13, 14]. As AMB is a non-linear system, other methods like feedback linearisation or sliding mode 
control can also be applied [10, 15]. With an effort to eliminate uncertainties in plant modelling, Bonffito et al 
propose an offset free control for AMB based on classical model predictive control [16]. 
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It is evident that in previously mentioned research, hard limitation on control input and output is 
omitted [17-20]. Control saturation might result in control degradation and violation of output constraint leads 
to system mechanical failure. In this paper, we have adopted backstepping control algorithm to regulate and 
stabilize the operation of AMB. The backstepping method há been employed in robotics [21], process control 
[22], space applications [23, 24], and in AMB systems [25]. It is proven to be suitable with strict-feedback 
system and to have the flexibility of removing instability while avoiding cancellation of potentially useful 
nonlinearities [21]. The contribution of the paper can be named is the consideration of bounded system input 
and output in control design. This paper is organized as the following. The mathematical model of AMB is 
first developed in Section 2. Controller design process is presented in Section 3. In Section 4. simulation results 
are provided together with the discussions. Finally, Section 5. concludes this paper. 


2. MATHEMATICAL MODEL OF A DUEL COILS MAGNETIC ACTUATOR 

It is fundamental that the magnetic force is proportional to the current square. Thus, regulating the 
current can result in the force change. It is assumed that the rotor shaft has already been elevated along z axis 
in vertical direction by another system. The system includes 2 pairs of the same electromagnets along x and y 
axes in horizontal directions, ones of each pair are placed in the opposite position as illustrated in Figure 1. 
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Figure 1. Four electromagnet system 


Each pair, then produce the forces of attraction F&F and F3&F which are adjusted by regulating 
the currents 21, 22,23 and 24 respectively so that the shaft can be kept balance in the space within those magnets. 
Assuming that (x1, 71), (£2, 72), (£3, i3), (£4, i4) are the positions and currents of electromagnets 1, 2, 3, and 4 
respectively. The expression of the magnetic forces is given as 


: : 2 z ; 2 
Fi = HN?’ Ag = 5 (5) , Fo = Higa Ag 7 7 (2) 


4g? 4 Ly Ars? 4 D2 (1) 
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where K is a coefficient and calculated as K = gN? Ag, where ug is the permeability of air, N is 

the number of turns in each coil, and A, is the cross-section area of the electromagnet. It is assumed that the 

inertia and geometric rotating axes of the rigid rotor coincide to each other, hence, the central point G is the 

mass center of the rotor and m as its mass. its mass. The x axis direction forces (1) exerted on rotor result in 

translational and rotational motions such that x and 0, DOF’s force and torque equations are given as (2) and 
(3) respectively: 

m(Zg) = (Fy — F2) + (F3 — Fa) (2) 
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1.61.00, + (Fi — Fo).Du — (F; — F4).Di (3) 


Where J, and I, are the total moments of inertia about axial and radial direction axes z, x and y 
through the rotor’s mass centre or center of weigth, respectively, 0y is the rotor angle over y axis, D,, is the 
distance form the top electromagnets to rotor central G, D; is the distance form the bottom electromagnets 
to rotor central G, —J tO is the reinforced torque of rotor, as shown in Figure 1. If the distance between 
rotor and the magnet at stable position is x9, movement of rotor within the top two magnets is x,, and that 
within the bottom two magnets is xı, the distances £1, £2, %3, x4 between rotor and each magnet can be calcu- 
Lı = Qo — Ly 


and with respect to the bottom two magnets: 
9 bs to 0 ee aed ON 


lated as with respect to the top two magnets: 


3 = Xo =T l l ) 
l i a 4 7 . It is found that the AMB 4th order system with two pairs of electromagnets arranged as in 
4 = To F Tı 


Figure 1 can be separated into two magnet systems, or simplified to two 2nd order systems. In that case, with 
an assumption of very small 6,,, the movement of rotor is represented as, wrt. the two upper magnets: 


Uy = fo + De, (4) 


and for the two lower magnets: 
Ll = Lg — D6, (5) 


Taking 2nd order derivatives of (4), and combining with (2) and (3), we have 


Ëu = Ëg + Diy 


1 1 1 1 
== fy pf Dy (AD. ERD.) 
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m I, L; (6) 
a aoe E il 
és (£o — Zu)? (£o + Ly)? 
where a, = Bu (4 + au). It is noted that the coupling term related to 6, and 6, is omitted. 
The coupling effects is considered as system disturbances. Applying the same procedure, from (5): 
%, = 4, — D6, 
2 2 (7) 
| rn aa 
(zo = 21) (Xo + £1) 
where a; = ae (4 + PE). On the other hand, applying Kirchhoff’s voltage law for each coil, we 
have the following equations: 
l di, Kd/ù , dig Kd (‘2 
Se a i a g 
ee e+ea(s) Sone aia (2) (8) 
, di Kd [is diy) Kd (% 
=R Lig Tae aa Nt oa, AS =R L, — — — | — 9 
a eas TERT ee e+ea(=) ”) 
Deriving form (8) and (9), the currents are represented as: 
2(xo — Lu) K.vyt1 2(£o + Lu) K.vyte 
h= (u — Ray — ———— ), = (tn — Rig + ——— 5), (10 
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In (4), (10), together presents the mathematical model of the two upper magnets as below: 


Li Üy 
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Similarly, (5), (11) provides the mathematical model of the two lower magnets: 


£i = V] 
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In summary, the system model of AMB studied in this work consists of (12) and (13) and will be used 
in the subsequent sections to design controller and investigate the operation. 


3. CONTROLLER DESIGN 
It is assumed that the speed can be estimated as derivative of position, the effect of rotor speed on 
system operation is negligible, and magnetizing currents are taken as control input. 


3.1. Control law of the two upper electromagnets 

Step l: Find the controller to enable the position x,, of the rotor track desired set point at the stable 
value which is O (along Ox axis). If zı is the difference between the rotor position and the stable one: z1 = £u, 
and its derivative is 2] = £u = Vu. Considering the following barrier Lyapunov candidate function: 


ks 


aie 14 
y (14) 


1 
Vi= 5m 


where ky is the limit of z1. It is clear that V; (z1) is radially unbounded as zı approaches ky or —kp. 
The Barrier Lyapunov function is used to reduce the error in rotor shaft position when comparing with the 
desired value, so as it would prevent the rotor shaft move too far away with a large distance which is greater 
than the air gap. This would lead to the collision between the rotor and the magnets, then damage the system. 
The derivative of (14) is 
B 21-2] 
ke 2 





Vj (15) 


Based on Lyapunov stability, it is required that V, < 0, thus virtual control function can be selected as 


Vudk = — (kg — zi) kızı (16) 
where kı is a positive constant. Then, V; = TrA = —kız? < 0 satisfies the stable condition. 
b 1 
Let vuak = 1, we have: 
: ‘ Oa ‘ 2 2\ à 
Uidie =O = ayo = (—ky kj a 3k 27) ZY (17) 
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Thus, va is the virtual control which facilitates x,, reach the set points. 
Step 2: Identify virtual control to regulate v,, to match vua. If deviation of vu from Vudk 1S 22: 


Lo. = Uy. ye = i = On (18) 


Or it can be represented as v, = Vudk + Z2. Derivative of (18) results in: 


Oa, 


=e jl 
a (19) 


The Lyapunov candidae function in this step is chosen as: V2 = V; + a: we then differentiate both 
side to get: 


: 2 Z122 . ; 
Vo = —kı zí + o TO= a1) (20) 
b 1 
In order to have i> < 0, the virtual control function is selected as ùudk = Q2 = — kə Z2 + &1 — mee 
b 1 


where kz is a positive constant. Substitute ùuak in (3.1.) for ù, in (20) , we have: 


ZIZ . va ; 
P 5 + 22(—kaze + Gy — s —å1) 
kg — 24 ky — 23 (21) 


2 2 
= —ky 2zi — k225 


V = -kız7 + 


In (21) shows that Vo < 0 as required for stability. Therefore, v,, as virtual control law is identified. 
It is a function of 2; and tə based on (12): 








; 2 ; 2 
i te ee N, (22) 
LO Ür Lect ey 
Let ay, = Q2/dy, from the above equation, it can be shown that a, = Pou zy + a Un 


Step 3: Design the current control law such that current 7 would match the set point tą. As presented 
in Step 2, the virtual control law ù, is a function of 2 currents 21 and 22, which are equivalent to electromagnetic 
forces of the two magnets. The fact that these two magnets operate simultaneously to maintain electromagnetic 
forces leads to higher energy consumption. Thus, a control scheme of switching on and off the two currents 
sequentially is employed to achieve energy savings as the following: 

Case 1: x, < 0 and 22 = O, it is shown that 


iid = (£o — Lu) VQu (23) 











Oz1 Ov Lo— Tu 
between 2; and set point tiq, 1.€.: Zy1 = îi — g Differentiating both side of Equation 3.1., we get: 2,1 = 


i1 — 11. In order to limit the input signal, i.e. current, within a bounded range, the current variable is provided 


aS 
ome 24 
AER (=) (e 


where im is the magnitude of current range, and v is the coefficient of tanh(). Consider the Lyapunov 
candidate function in this step as: V3 = V2 + ane the derivative of this Equation is: 


; ! ! ; 2 
On the other hand: ijg = 242, + 244%,,, where ù, = ay. ( al ) . Call 2,1 is the deviation 








Ji Ji 
Va eee hks s (25) 
Oz OVa 
Based on (25) and the condition that Va < 0, the control law J, is selected as 
Oiid. Oiid. 
Ty = -kyo +t Eia (26) 


Z1 
oz OVy 
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where k,,; is a positive constant. Substitute this /; in (25), we have 


: , Oiid. Olid . Olid. Oiid. 
V3 =V: vl | —hy1Zy a, ee: i oa u 
3 at en ( ivi Tt Ta eE ) z ži Di. (27) 








_ T 2 
=V kyz 


It can be seen from (27) that V3 < 0 which satisfies stable condition. Thus, with x,, > 0, /; as in (26) 
is the control law to stabilize the upper part of the rotor. 
Case 2: x, > 0 with respect to i; = 0. The condition implies that: 22g = (£o + £u) V—&u where 





2 : 
Urm nE ee (=) . Call z,2 is the deviation between tə and set point 22g: Zv2 = 12 — i2q. Differentiating 


both side of (3.1.), we get Z,2 = i= isq. Similar to case 1: 


ig = Ip r 


The Lyapunov candidate function in this case is V4 = V2 + ee The derivative of this Lyapunov 
function is: 








. Ji Ji 
Vea dee ue an (29) 
Oz OVa 
With the condition of Va < 0, the control function [5 is selected as: 
Oid . Oleg . 
Ip = -kizo + = 41 + qty (30) 
Oz OV, 
where k,,2 is a positive constant. Substitute the selected /> in (29), we have: 
: i Oizd . Olod . Olod . Ol2d . 
V4 =V + 2y2 —ky2 22 TE A F i z A = ai 
Oz OVy Oz OVy, (31) 


_ T 2 
= Vo = hizo 


In (31) shows obviously that V; < 0 which satisfies stable condition, and the control law I» selected 
can stabilize the upper part of the rotor. 


3.2. Control law of the two lower electromagnets 

The design procedure is similar to that of two upper electromagnets as presented in a). It also includes 
3 steps as the following: 

Step 1: Identify position control xı to reach the stable position, which is O (along Ox axis). Let z3 be 
the deviation between rotor shaft and the stable position, 1.e.: 23 = 41 > 23 = ťčı = vı The barrier Lyapunov 


: . : k2 : f . 
candidate function is V5 = 1 ln po The virtual control is chosen as Vidak = — (k? — z3) k3z3 where k3 is a 
b 3 


positive constant. Similarly, it can be proven that this control law renders V; < 0. Let vigx = a3, and compute 
its time derivative 


: Oa , 
Vldk = Q3 = Jz = (—k3ky + 3k325) Z3 (32) 


Step 2: Select virtual control so that v; would be able to reach viag. Let the difference between v; and 
Vidk De Z4: Z4 = VU — Vidk = Vi — a3. Or Y = Viak + 23 In this step, the Lyapunov candidate function is 
Ve = V; + oe We pick the virtual control to satisfy that Vg < 0 as: 


, z 
Didk = 04 = —kaza + G3 — = 5 (33) 
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where k is a positive constant. In (13) provides the calculation of ù; from 23 and 74 as: 


; 2 : 2 
D =a( 2 ) -a = ) (34) 
£o — XI £o + Tı 


Let a; = &3/aı, we have a; = Fal ag + Sd 


Step 3: The switching scheme of currents supplied to lower magnets are 


Case 1: x; < 0 and i4 = 0 implies that: i3g = (£o — x1) ,/a; And thus, isq = Fa Zg + sa in, Let 2,3 is the 


deviation of 23 from set point 734, we have: 2,3 = t3 — 73q. Using tanh to limit 23 in the required range: 











is = Ís 
ig = im tanh (2) (35) 


tm 


Barrier Lyapunov function in this step is V7 = Vg + a In order to render V- < 0, virtual control I3 is 
selected as 


Oia, . 
p E se (36) 
Oz3 Ov, 


where k,,3 is a positive constant. 
Case 2: x; > 0 and 73 = 0 yields i44 = (£o + £1) V—a; and its derivatives is i4g = Gad Zg+ Sha Ùl, 











, 2 
where ù, = —a). (i ) . Applying the analogous design. then selecting virtual control 74 is shown as below 
lig Otad . 
Iy = —hyazya + S23 + S01 (37) 
Oz3 Ov, 


where k,„4 is a positive constant. 


4. SIMULATION AND DISCUSSION 

Numerical simulation parameters used in the study are presented as: Rotor mass m=5kg; number of 
coil turns V=400 turns; nominal air gap 79=0.001m; maximum position error k,=0.001m; initial position of 
upper rotor shaft x,,=0.0001m; initial position of lower rotor shaft x;=0.0001m; self inductance L, is 0.001H; 
cross section area of iron core A, is 0.001m?; permeability of air gap j1,=1.256 x 10-°H/m; moment of inertia 
I,=2.900 x 10~?kg - m?; distance from rotor central to upper magnets D,,=4.166 x 10~*m; distance from 
rotor central to lower magnets D,=7.602 x 10~2m. Controller’s coefficients are of kı = 11; kə = 1700; ky, = 
700; k2 = 10000; k3 = 10;k4 = 1600;k,3 = 700;k,4 = 10000. In the paper, to emphasize the ability of 
handling input and output constraints of the proposed controller, the rotor shaft is driven to equilibrium position 
and accelerating to 1000rpm. This simulation procedure implies the effects of coupling term related to 6, and 
0y can be eliminated. 

Case study 1: current limit is îm = 3A. As shown in Figure 2a and 2b, the upper and lower body of 
the AMB can be regulated from its deviation to the stable position within 0.01 second. Duration to reach the 
zero displacement lower body is also around 0.01 second, however the overshoot is a little bit more, 1.e. around 
5pm, it is clear that the value is well below the threshold define by kẹ. Meanwhile the duration of central 
displacement is corrected within the same interval and the overshoot is slightly smaller than that of lower body. 

Case study 2: current limit is îm = 2A. The current limit is reduced to 2A in this case, but the initial 
displacement of the rotor shaft is kept the same. It is clearly observed that the settling times of the upper 
body and lower body in Figure 2a and Figure 2b are slightly longer that those in case one in Figure 3a, 3b 
respectively. It 1s due to the fact that controllers need to take more effort to stabilise the system with smaller 
current fed thanks to the use of the hyperbolic tangent function in the design. The peak currents supplied to the 
AMB electromagnetics are all less than the provided limit as shown in both two cases as illustrated in Figure 4 
and Figure 5. These peak values for lower magnets are also less than those of the upper ones. When the current 
limit is decreased, it is observed that there is more oscillation of current response. The cause can be explained 
as less magnetic forces provided to the system. 
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Figure 2. Rotor displacement, (a) Upper body displacement, and (b) Lower body displacement 
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Figure 3. Rotor displacement, (a) Upper body displacement, and (b) Lower body displacement 
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Figure 4. Current responses, (a) Upper magnets, and (b) Lower magnets 
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Figure 5. Current responses, (a) Upper magnets, and (b) Lower magnets 
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5. CONCLUSIONS 

In this paper, a 4th order AMB has been modeled as two 2nd order subsystem with magnetizing cur- 
rent is treated as control input. The backstepping method is adopted in control design for the obtained model. 
The controllers have been built and validated via simulation in different case studies in a view of input satu- 
ration and bounded output . It is shown that our proposed approach is able to facilitate the AMB regulate gap 
deviations as desire and thus stabilizes the system. Future work include practical implementation of the whole 
system, it would enable further investigation of the proposed works thoroughly for real-life applications. 
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